home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / MATHEMAT / 1769.ZIP / ALGIN3 < prev    next >
Text File  |  1987-02-21  |  730b  |  37 lines
  1. order e p q pb qb
  2. *  e is a small parameter (epsilon)
  3. define h
  4.   p^2
  5.   q^2
  6.   8 e p q
  7.  
  8. define f
  9.   pb q
  10.   e pb^2
  11.   -e q^2
  12.  
  13. differentiate dfdq = d f/ d q
  14. differentiate dfdpb = d f/ d pb
  15. print dfdq
  16. print dfdpb
  17. *  There is no "solve" command in PFSA, one can use various
  18. *  commands such as "compute" and "extract" to aid in big
  19. *  computations.  Here we see by a simple calculation:
  20. *  q=qb - 2 e pb  and p=pb - 2 e qb + 4 e^2 pb .
  21. define qequals
  22.   qb
  23.   -2 e pb
  24.  
  25. define pequals
  26.   pb
  27.   -2 e qb
  28.   4 e^2 pb
  29.  
  30. substitute qequals for q in h
  31. substitute pequals for p in h
  32. print h
  33.  
  34. *  We see that the e term was eliminated.  More transformations
  35. *  can be done to eliminate other unwanted terms.
  36. end input
  37.